Draft version February 4, 2008 

Preprint typeset using I4TgX style emulateapj v. 03/07/07 



THE FINAL SPIN FROM THE COALESCENCE OF ALIGNED-SPIN BLACK-HOLE BINARIES 
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ABSTRACT 

Determining the final spin of a black-hole (BH) binary is a question of key importance in astrophysics. 
Modelling this quantity in general is made difficult by the fact that it depends on the 7-dimensional space 
of parameters characterizing the two initial black holes. However, in special cases, when symmetries can be 
exploited, the description can become simpler. For black-hole binaries with unequal masses but with equal 
spins which are aligned with the orbital angular momentum, we show that the use of recent simulations and 
basic but exact constraints derived from the extreme mass-ratio limit allow to model this quantity with a simple 
analytic expression. Despite the simple dependence, the expression models very accurately all of the available 
estimates, with errors of a couple of percent at most. We also discuss how to use the fit to predict when a 
Schwarzschild BH is produced by the merger of two spinning BHs, when the total angular momentum of the 
spacetime "flips" sign, or under what conditions the final BH is "spun-up" by the merger. Finally, suggest 
an extension of the fit to include unequal-spin binaries, thus potentially providing a complete description of 
the final spin from the coalescence of generic black-hole binaries with spins aligned to the orbital angular 
momentum. 

Subject headings: black hole physics - relativity - gravitational waves - stars: statistics 



1. INTRODUCTION 

The determination of the final spin of a BH binary is a 
question of key importance in astrophysics. Modelling this 
in general is made difficult by the fact that it depends on the 
7-dimensional space of parameters characterizing the two ini- 
tial BHs. However, in special cases, when symmetries can be 
exploited, the description can be much simpler. 

Several recent studies have shed light on the remnant of the 
merger pro cess. Using conservation pr inciples, Hughes and 
Blandford dHughes & Blandfordll2003|) argued that mergers 
rarely lead to rapidly rotating objects. IGonzalez et al.l ( l2007al) 
numerically evolved a sequence of non-spinning unequal- 
mass BHs, arriving at detailed estimates of the r adiated energy 
and angular momentum. In a series of papers dKoppitz et al.l 
120071: iPollnev et al.ll2007l:lRezzolla etai]|2007l) we have stud- 
ied the parameter space of mergers of equal-mass BH bina- 
ries whose spins are aligned with the orbital angular momen- 
tum but otherwise arbitrary. The findings agree we l l with 
independent numeric al evolutions (ICampanelli et al.l 120071 : 
iHerrmann et al.ll2~007l) . as well as more recen t studies of mod- 
els w ith initial spins up to J/M 2 = 0.8 (M arronetti et al.l 
|2007|) . An important result of these studies has been the deter- 
mination of simple (quadratic) fitting formulas for the recoil 
velocity and spin of t he merger remnant as a function of the 
initial BH parameters dRezzolla et al.ll2007l) . 

A number of analytical approaches have been developed 
over the yea rs to determine the final spin of a binary 
coalescence dPamourl 1200 U IBuonanno & Damourl |2000|: 



Buon anno. et al. 1 120061: iDamour & Nagarl 120071: iBovle etail 
20071) . Very recently, an interesting method, inspired by the 



dynamics of a test particle a round a Kerr BH , has be en pro- 
posed for generic binaries (IBuonanno et al.l (l2007bl) . BKL 

1 Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, 
Potsdam-Golm, Germany 

2 Center for Computation & Technology, Louisiana State University, Ba- 
ton Rouge, LA, USA 

3 Department of Physics and Astronomy, Louisiana State University, Ba- 
ton Rouge, LA, USA 



hereafter). The approach assumes that the angular momen- 
tum of the final BH is the sum of the individual spins and of 
the orbital angular momentum of a test particle on the last- 
stable orbit of a Kerr BH with the same spin parameter as that 
of the final BH. 

Here, we combine the data obtained in recent simula- 
tions to provide a phenomenological but analytic estimate 
for the final spin in a binary BH system with arbitrary 
mass ratio and spin ratio, but in which the spins are con- 
strained to be parallel to the orbital angular momentum. 
Our numerical s imulations have bee n carried out using the 
CCATIE code dPollnev et al.l \2001b. In addition to the 
data presented in iRezzolla et al.l d2007l) , we add three sim- 
ulations of equal-mass, high-spin binaries and three sim- 
ulations of unequal-mass, spinning binaries (see Table [T]). 
Oth er data is taken from unequal-mass, nonspinning bina- 
ries dGonzalez et al]l2007at iBerti et aTJl2007t Buonann o et alJ 
2007al). and of equal-mass, sp inning binaries dRezzolla et alJ 
20071: IMarronetti et al.l 120071) : all of the AEI data is sum- 
marized in Table Q] To avoid the possible contamination 
fro m the errors associated with high-spin binaries reported 
by IMarronetti et alJ d2007l) . we have not considered bina- 
ries w ith initial spin | J/M 2 \ ^ 0.75 reported in the liter- 
ature (ICampanelli et alJ l2007t IMarronetti et al.1 120071) . We 
have, however, considered estimates of high-spin binaries 
(cf., Table [TJ, for which we know the spins remain es- 
sential ly constant prior to merger, with changes less than 
0.5% dPoHney et al.ll2007l) . and that are very well captured 
by the fit. 

2. METHODS AND RESULTS 

We start by considering the final spin Ofi n as a function of 
the two free variables in the problem: the symmetric mass 
ratio v = M\M 2 j(M\ + M 2 ) 2 and the spin of the initial BHs 
a = J/M 2 , i.e., af\ n = Jfi n /M| n = afin(a, v). (Note a is 
dimensionless and not the angular momentum per unit mass.) 

By construction a\ = a 2 = a, and a/\a\ = ±L/\L\, where 
L is the orbital angular momentum. We next express ag n as a 
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TABLE 1 

Initial parameters of the new binaries computed at the AEI. 
The different columns contain the initial spin a, the 
symmetric mass RATIO V, half of the initial separation 

x/M = i(xi — X 2 ), THE DIMENSIONLESS INITIAL ANGULAR 

MOMENTUM J = J/(flM), THE NUMERICAL AND FITTED VALUES FOR 
a fln AND THE CORRESPONDING RELATIVE ERROR. 





a 


V 


x/M 


J 


a ti„ 


fin 


|err.|(%) 


t8 


-0.5840 


0.2500 


3.1712 


2.432 


0.4955 


0.4981 


0.53 


ta8 


-0.3000 


0.2500 


3.7078 


3.000 


0.5941 


0.5927 


0.23 


m 


-0.8000 


0.2500 


3.8082 


2.200 


0.4224 


0.4227 


0.08 


mi 


-0.8000 


0.2500 


4.8600 


2.400 


0.4266 


0.4227 


0.92 


P i 


-0.8000 


0.1580 


3.2733 


0.336 


0.0050 


0.0046 


9.89 


p-2 


-0.5330 


0.1875 


3.3606 


1.872 


0.2778 


0.2794 


0.57 


P 3 


-0.2667 


0.2222 


3.4835 


2.883 


0.5228 


0.5216 


0.23 




third-order polynomial of v and a 

ctfin = So + s i a + s 2a 2 + s 3 a 3 + s±a 2 i> + s^av 2 + 

t au + t x v + t 2 v 2 + t 3 v 3 . (1) 

Expression (HJ is a lowest-order ansatz. It intends to capture 
the behaviour of a function known exactly only in the extreme 
mass-ratio limit (EMRL) and which has support from numer- 
ical simulations in two restricted regimes: i.e., v = 1/4; < 
\a\ < 0.75 and 0.16 < v < 1/4; a = 0. A-priori there is 
no reason to believe expectation that afi n (^, a) is that the pro- 
posed fit will capture the general behaviour well, but in fact it 
does. 

Given the available numerical estimates, it is possible to 
calculate the coefficients so-s.5, and to-t^ by simply perform- 
ing a two-dimensional (2D) least-square fit of the data. This, 
however, would require a lot of care and is likely to lead to 
inaccurate estimates. This is because the space of parameters 
presently accessible to numerical simulations is rather small. 
Reliable results are in fact available only for spins \a\ < 0.8 
and mass ratios q = M2/M1 > 0.25 and thus corresponding 
to v > 0.16. However, it is possible to exploit exact results 
which hold in the EMRL, i.e., for v = 0, to constrain the 
coefficients in expression ([TJ. It is worth emphasizing that 
the EMRL results are not only exact, but also in regimes that 
numerical relativity simulations cannot probe. More specifi- 
cally, we can exploit that in the EMRL the final spin cannot 
be affected by the infinitesimally small BH. In practice, this 
amounts to requiring that 

Ofi n (o, v = 0) = a , (2) 

which constrains four of the six coefficients 

s = s 2 = s 3 = , si = 1 . (3) 

Additional but non-exact constraints can also be applied by 
exploiting the knowledge, near the EMRL, of the functional 
dependence of Ofi n on the mass ratio. A convenient way of 
doing this is suggested by BKL, and within this approach we 
perform a Taylor expansion of <Zfi n for v -C 1 and determine 
that 

4ml (0=1,^=0) = 20/3/3 - 1) , agj (o=0il/=0) = 2^3 , 
<n\ {a =-l,u=0) = 2 (1 + 19>/l5/45) , (4) 



FIG. 1. — Global dependence of the final spin on the symmetric mass 
ratio and on the initial spins as predicted by expression {3). Squares refer to 
numerical estimates while circles to the EMRL constraints. 

where a' Rn = da^ n /dv. The coefficients in (Q]i are 
then s 4 = V3(19\/5 - 75)/45, h = 2^3, t = 
[V3(15 - 19^) - 90] /45. While this may seem a good 
idea, it leads to bad fits of the data. We believe this is due 
to two distinct reasons: (i) the lack of accurate numerical 
data for near-extreme BHs, i.e., \a\ « 1, and which there- 
fore leads to incorrect estimates of the coefficients; (ii) ex- 
pressions (HJi are analytic but not exact and should be used 
with caution. There are, in fact, deviations from analyticity in 
v as v — > 0, and as revealed by the presence of integer powers 
of v 1 ' 5 during the transition between t he las t stable orbit and 
the plunge (see iBuonanno & Damourl d2000l) ). In the case of 
non-spinning binaries (a = 0), it is n ow possible to verify tha t 
the deviations are indeed very small (iDamour & Nagarl2007l) . 
but this check is not possible for very large spins. In view of 
this and to make the minimal number of assumptions, we re- 
tain the analytic estimate only for the coefficient t\, so that (fl~|i 
has five out of ten coefficients constrained analytically 

flfin = a + s±a 2 v + s^av 2 + t$av + 2\fiv + t2V 2 . (5) 

Determining the remaining five coefficients from a least- 
square fit of the available data yields 

s 4 = -0.129 ± 0.012 , s 5 = -0.384 ± 0.261 , 
t = -2.686 ± 0.065 , t 2 = -3.454 ± 0.132 , 
t 3 = 2.353 ± 0.548 , (6) 

with surprisingly small residuals and large error-bars only for 
S5. The functional behaviour of expression (|5]l and the posi- 
tion of the numerical data points are shown in Fig.Q] 

In the following we discuss the properties of the proposed 
fit, providing evidence that it represents a very accurate de- 
scription of the available estimates, and discuss how to use it 
to make astrophysically interesting predictions. 

( i) The estimate for the final spin in the case of equal masses 
and the comparison with available data and estimates is made 
in Fig. [2] The upper pa nel shows the numeri cal estimates, 
[circles for the A EI data (Rezzolla et al. 2007) and stars for 
the FAU-Jena data dMarronetti et alj2007l) l. the BKL estimate 
and our 2D fit through ©. The lower panel shows the resid- 
uals between the different estimates and the 2D fit; these are 
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FlG. 2. — Upper panel: Comparison of the numerical data with the 2D 
fit through (3} in the case of equal-mass binaries, (v = 1/4). Empty 
circl es indicate the AEI d ata (Rezzolla et al. 2 0071) . stars the FAU-Jena 
data ( MaiTonetti et al. 2007)], a long-dashed line the BKL, and a short-dashed 
one the fit. Lower panel: residuals between the different estimates and the fit. 
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FlG. 3. — Upper panel: Comparison of the numerical data with the 2D 
fit throug h (5) in the case o f nonspinning binaries. Empty circles indicate th e 
Jena data I Berti et al. 2007), stars the Go ddard data I Buonanno et al. 2007af) l, 
a long-dashed line the quadratic EOB fit I Damour & Nagar 2007) and a short- 
dashed line our 2D fit. Lower panel: residuals between the different estimates 
and the 2D fit. 



always of a few percent only and become larger for the BKL 
estimate when a < 0. 

(ii) Despite the cubic dependence assumed in (Q~|), expres- 
sion (0 is only quadratic wi th a. When v = 1/4 , it confirms 
what was obtained recently dRezzolla et alJl2007). indicating 
that, for equal-mass binaries, the next order will be four. 

(Hi) Using ([5]) and (0 we estimate that the minimum and 
maximum final spins for an equal-mass binary are aa n = 
0.3502 ± 0.03 and a fi „ = 0.9590 ± 0.03, respectively. 

(iv) For nonspinning binaries, expression (0 is cubic in v 
and a comparison with the available data and the estimate 
from the EOB approach c ombined with test-mass limit pre- 
dictions for the ringdown dDamour & Nagar 2007) is shown 
in Fig. [3] In particular, the upper panel shows the nume rical 
values, [empty circles for th e Jena data dBerti et al.l 20071) and 
stars for the Goddard data dBuonanno et alJl2007a )1, a long- 
dashe d line for the quadratic EOB ID fit dDamour & Nagarl 
120071) and a short-dashed line for our 2D fit. (Because it is 
very similar to the EOB estimate, we have not shown the BKL 
prediction.) The residuals in are shown in the lower panel. 

(v) A physically useful condition that can be deduced from 
the 2D fit are the values of the initial sp in and mass ratio that 
will lead to a final Schwarzsc hild BH dHughes & Blandfordl 
l2003t iBuonanno et al.l l2007bl) . In practice this amounts to 
requiring afi n (a, v) = in © and this curve in the (a, v) 
plane is shown in the upper panel of Fig. |U Binaries on the 
curve produce Schwarzschild BHs, while binaries above the 
curve start with a positive total angular momentum and end 
with a positive one; binaries below the curve, on the other 
hand, start with a positive total angular momentum and end 
with a negative one, i.e., with a global flip. Also shown 
in the upper panel of Fig. [4] is the prediction from BKL: 
a schw Ibkl = 2uy/3/(2u — 1). The two estimates are very 
similar for all values of v and small differences appear for 
v > 0.15, where the BKL estimate is less accurate. Shown 
with a cross is the binary p\ (cf., Table[T|i which yields a final 



BH with spin a^ n = 0.005. The numerical value is between 
the BKL prediction and the 2D fit. 

(vi) The BKL is expected to be particularly accurate for 
v <C 1 and its prediction in this regime are captured very well 
by the 2D fit (of course the two predictions are identical for 
v = 0). This is shown in the lower panel of Fig.@]with differ- 
ent curves referring to v = 0.001 , 0.01 and 0.1; interestingly, 
the differences are small even for v = 0.1. 

( vii) It is simple to derive the value of a which will produce 
a final BH with the same spin as the initial ones. This amounts 
to requiring that afi n (a, v) — a in (0 and the resulting solu- 
tion is shown in Fig. |5J clearly, the axis v = is a trivial 
solution and a magnification of the behaviour away from the 
EMRL is shown in the inset. For equal-mass binaries the crit- 
ical value is a CI j t = 0.9460, in very good agreement w ith the 
BKL estimate a crit > 0.948 (IBuonanno et all l2007bl) . The 
minuteness of the region for which afi n < a (dashed region) 
suggests that BHs from aligned-spins binaries are typically 
spun-up by mergers. 

(viii) It is easy to verify that by setting v = 1/4 and 
2a = ai + ci2 in ©, the coefficients S1-S5 and £n-*3 coin- 
cide, within t he error-bars, with th e coefficients po, p\ and 
P2 reported in lRezzolla et al.l (120071) for equal-mass, unequal- 
spin binaries. The fact that the fit here is equivalent to, 
but has been independently derived from, the one for the 
equal-mass, unequal-spin binaries, is an indication of its ro- 
bustness. Indeed, it is possible to extend (0 to the whole 
(01,02,1^) space i.e., to describe the final spin of generic 
aligned, unequal-spin, unequal-mass BH binaries, by replac- 
ing a with (ai + a2<7 2 )/(l + q 2 ). The resulting expression re- 
duces to (T5j> for unequal-ma ss, equal-spin binaries, and to the 
one in lRezzolla et al.l d2007l) for equal-mass, unequal-spin bi- 
naries. Our suggested extension of © to the (01, 02, v) space 
is the simplest one which recovers, for aligned spins, the well- 
tested limits of equal-mass, unequal-spins and unequal-mass, 
equal-spins. Work is in progress to validate this ansatz with 
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A final comment is one of caution. The dependence of the 
final spin on the mass ratio in the case of extreme aligned 
BHs is particularly challenging to calculate and not yet inves- 
tigated accurately by numerical calculations. The predictions 
of expression (0 in this limit amount to mere extrapolations 
and are therefore accurate to a few percent at most. As an 
example, when a = 1, the fit (0 is a non-monotonic function 
with maximum a^ n ~ 1.029 for v ~ 0.093; this clearly is an 
artifact of the extrapolation. 
-1 -0.8 -0.6 = -0.4 -0.2 
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FIG. 4. — Upper panel: Set of initial spins and mass ratios leading to a 
final Schwarzschild BH: i.e., afl n (a, v) = 0. The two curves refer to the 
BKL estimate (long dashed) and to the 2D fit (short dashed), respectively. 
Indicated with a star is a numerical example leading to aa n = 0.005. Lower 
panel: Comparison between the BKL prediction (symbols) and the 2D fit 
(solid, dashed and long-dashed lines) near the EMRL. Different curves refer 
to different values of v and the match is complete for v = 0. 
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FIG. 5. — Critical values of the initial spin and mass ratio leading to a final 
BH having the same spin as the initial ones i.e., afi n (a, v) = a. A magni- 
fication is shown in the inset, where the dashed/non-dashed region refers to 
binaries spun-down/up by the merger. 



3. CONCLUSIONS 

Modelling the final spin in a generic binary BH merger is 
not trivial given the large space of parameters on which this 
quantity depends. We have shown that the results of recent 
simulations combined with basic but exact considerations de- 
rived from the EMRL allow us to model this quantity with a 
simple analytic expression in the case of BH binaries having 
unequal masses and unequal spins which are aligned with the 
orbital angular momentum. When compared with all other 
estimates coming either from numerical calculations or from 
approximation techniques, the estimates of the 2D fit show 
differences which are of few percent at most. 
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